Symmetric BEM formulations for 2-D steady-state and transit potential problems

This paper presents a new concept for symmetric boundary element method (SBEM) applicable to 2-D steady-state and transit potential problems. Two kinds of SBEM formulations are derived. Symmetry is obtained simply through matrix manipulation, and no hypersingularity appears. Therefore, SBEM is much easier than the traditional symmetric Galerkin BEM. Compared with the traditional asymmetric BEM, the present SBEM can reduce the computational cost for time domain problems only. However, when applied to BEM/FEM coupling procedure, SBEM can reduce the computational cost for both steady-state and time domain problems. Three numerical examples are included to illustrate the effectiveness and accuracy of the present formulations.     

Transition to turbulence in rotating-disk boundary layers—convective and absolute instabilities

This work is an experimental study of mechanisms for transition to turbulence in the boundary layer on a rotating disk. In one case, the focus was on a triad resonance between pairs of traveling cross-flow modes and a stationary cross-flow mode. The other was on the temporal growth of traveling modes through a linear absolute instability mechanism first discovered by Lingwood (1995, J Fluid Mech 314:373–405). Both research directions made use of methods for introducing controlled initial disturbances. One used a distributed array of ink dots placed on the disk surface to enhance a narrow band of azimuthal and radial wave numbers of both stationary and traveling modes. The size of the dots was small so that the disturbances they produce were linear. Another approach introduced temporal disturbances by a short-duration air pulse from a hypodermic tube located above the disk and outside the boundary layer. Hot-wire sensors primarily sensitive to the azimuthal velocity component, were positioned at different spatial (r,θ) locations on the disk to document the growth of disturbances. Spatial correlation measurements were used with two simultaneous sensors to obtain wavenumber vectors. Cross-bicoherence was used to identify three-frequency phase locking. Ensemble averages conditioned on the air pulses revealed wave packets that evolved in time and space. The space–time evolution of the leading and trailing edges of the wave packets were followed past the critical radius for the absolute instability, r _{c A} . With documented linear amplitudes, the spreading of the disturbance wave packets did not continue to grow in time as r _{c A} was approached. Rather, the spreading of the trailing edge of the wave packet decelerated and asymptotically approached a constant. This result supports the linear DNS simulations of Davies and Carpenter (2003, J Fluid Mech 486:287–329) who concluded that the absolute instability mechanism does not result in a global mode, and that linear-disturbance wave packets are dominated by the convective instability. In contrast, wave-number matching between traveling cross-flow modes confirmed a triad resonance that lead to the growth of a low azimuthal number (n = 4) stationary mode. At transition, this mode had the largest amplitude. Signs of this mechanism can be found in past flow visualization of transition to turbulence in rotating disk flows.     

Non-linear current–voltage characteristics related to native defects in SrTiO3 and ZnO bicrystals

SrTiO₃ and ZnO bicrystals with various types of boundaries were fabricated in order to examine their current–voltage characteristics across single grain boundaries. Their grain boundary structures were also investigated by high-resolution transmission electron microscopy. In Nb-doped SrTiO₃, electron transport behaviors depend on the type of boundaries. Random type boundaries exhibit highly non-linear current–voltage characteristics, while low angle boundaries show a slight non-linearity. On the contrary, undoped ZnO does not exhibit non-linear current–voltage characteristics in any type of boundaries including random ones. It is suggested that the differences observed in current–voltage properties between the two systems are mainly due to the difference in the accumulation behavior of acceptor-like native defects at grain boundaries. A clear non-linearity is obtained by means of Co-doping even for the highly coherent Σ1 boundary in a ZnO bicrystal. This is considered to result from the production of acceptor-like native defects by Co-doping.     

电信运营商纵向一体化的经济学分析

阐述了我国电信服务业的纵向链条，以经济学的视角分析了电信运营商执行纵向一体化战略的动因．讨论了电信运营商纵向一体化的边界问题，最后提出我国电信运营商执行渐变一体化和战略联盟的建议。     

Exact integration in the boundary element method for two-dimensional elastostatic problems

This paper derives the exact integrations for the integrals in the boundary element analysis of two-dimensional elastostatics. For facilitation, the derivation is based on the simple forms of the fundamental functions by taking constant, discontinuous linear and discontinuous quadratic elements as examples. The efficiency and accuracy of the derived exact integrations are verified against five benchmark problems; the results indicate that the derived exact integrations significantly reduces the CPU time for forming the matrices of the boundary element analysis and solving the internal displacements.     

有限元数值试井解释技术在复杂断块油气藏评价中的应用

有限元数值模拟试井技术,可通过描述不同网格的瞬变压力对测试范围内各个小单元进行精细描述.综合运用解析试井、数值试井联带分析方法,将测试动态成果与地质静态成果有机结合、相互补充,可实现对复杂断块油气藏的精细描述,达到一口井搞清一个小型复杂断块油气藏的目的.     

An inverse solution methodology for estimating the diffusion coefficient of gases in Athabasca bitumen from pressure-decay data

An inverse solution methodology is developed for the estimation of diffusion coefficient of gases in highly viscous, oil-sands bitumens from isothermal, pressure-decay measurements. The approach involves modeling the rate of change in pressure using the diffusion equation for the liquid phase coupled with a mass balance equation for the gas phase. The inverse solution framework is utilized to arrive at two graphical techniques for estimating the diffusion coefficient. Both techniques involve the determination of the slope of a straight line resulting from plotting the experimental data in accordance with the developed model. An advantage of the proposed techniques is that the diffusion coefficient is estimated directly, i.e. without making it an adjustable parameter. The novelty of the proposed method is in its simplicity as well as its ability to isolate portions of the pressure-decay data affected by experimental fluctuations. The effect of the initial pressure on the predicted diffusion coefficient and pressure-decay profile was also investigated. The diffusion coefficients of CO₂, CH₄, C₂H₆ and N₂ in Athabasca bitumen at 50–90 °C and about 8 MPa were estimated and compared with literature values.     

地下水运动数值模拟过程中边界条件问题探讨

本文对地下水运动数值模拟过程中边界条件的涵义 和处理方法进行了分析和讨论。阐述了边界条件所包含的双重意义。指出随着人类活动影响 强度的日益增大，边界条件的处理要面临一些新的更为复杂的问题。在模型预报之前必须首 先对边界条件做出预报。边界条件的预报既要考虑自然因素的作用，同时也要考虑人类活动 (人工开采和人工补给)的影响及由于邻区水流条件变化而产生的耦合效应。之后，给出了两 个应用实例。     

Boundary knot method for 2D and 3D Helmholtz and convection–diffusion problems under complicated geometry

The boundary knot method (BKM) of very recent origin is an inherently meshless, integration‐free, boundary‐type, radial basis function collocation technique for the numerical discretization of general partial differential equation systems. Unlike the method of fundamental solutions, the use of non‐singular general solution in the BKM avoids the unnecessary requirement of constructing a controversial artificial boundary outside the physical domain. The purpose of this paper is to extend the BKM to solve 2D Helmholtz and convection–diffusion problems under rather complicated irregular geometry. The method is also first applied to 3D problems. Numerical experiments validate that the BKM can produce highly accurate solutions using a relatively small number of knots. For inhomogeneous cases, some inner knots are found necessary to guarantee accuracy and stability. The stability and convergence of the BKM are numerically illustrated and the completeness issue is also discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

A fictitious domain decomposition method for the solution of partially axisymmetric acoustic scattering problems. Part 2: Neumann boundary conditions

We present a fictitious domain decomposition method for the fast solution of acoustic scattering problems characterized by a partially axisymmetric sound‐hard scatterer. We apply this method to the solution of a mock‐up submarine problem, and highlight its computational advantages and intrinsic parallelism. A key component of our method is an original idea for addressing a Neumann boundary condition in the general framework of a fictitious domain method. This idea is applicable to many other linear partial differential equations besides the Helmholtz equation. Copyright © 2003 John Wiley & Sons, Ltd.